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Abstract 

Assuming that the maximal allowed number of identical particles 
in state is an integer parameter, q, we derive the statistical weight and 
analyze the associated equation which defines the statistical distribu- 
tion. The derived distribution covers Fermi-Dirac and Bose-Einstein 
ones in the particular cases q = 1 and q — > oo (rtj/q — > 1), respectively. 
We show that the derived statistical weight provides a natural combi- 
natorial interpretation of Haldane-Wu fractional exclusion statistics, 
and present exact solutions of the distribution equation. 
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1 Introduction 



Statistics which are different from Fermi-Dirac and Bose-Einstein ones be- 
come of much interest in various aspects. A recent example is given by 
Haldane-Wu fractional exclusion statistics (FES) [p], |2| which is used to de- 
scribe elementary excitations of a number of exactly solvable one-dimensional 
models of strongly correlated systems, and other models @, 0. This statistics 
is based on the statistical weight, which is a generalization of Yang- Yang |3J 
state counting as mentioned by Wu, 

_ (z i +(n i -l)(l-A))! 
1 n i \(z i -Xn i -(l-X))V {) 

where the parameter A varies from A = (Bose-Einstein) to A = 1 (Fermi- 
Dirac). This formula is a simple generalization and interpolation of Fermi 
and Bose statistical weights. While there is no physical meaning ascribed to 
A here, the physical interpretation of Eq. ([!]) is that the effective number of 
available single-particle states linearly depends on the number of particles, 

z/ = ^-(l-A)(n*-l), z\ = zi - X(m - 1), (2) 

for fermions and bosons, respectively. This is viewed as a defining feature of 
the fractional exclusion statistics. 

In the present paper, we show that the equation which defines Haldane- 
Wu statistical distribution can be derived from a different statistical weight, 
which has a clear combinatorial and physical treatment. Also, we present 
exact solutions of this equation. 



2 The combinatorics 

A number of quantum states of rii identical particles occupying Zi states, 
with up to q particles in state, 1 < q < rii, can be counted as follows. 

We consider a configuration defined as that it has a maximal possible 
number of totally occupied states (exactly q particles in state). A number 
of such totally occupied states is an integer part of rii/q which we denote by 
y . If q is a divisor of rii we have identically ^ = rii/ q, so that the number 
of unoccupied states is Zi — — . If q is not a divisor of rii we have one partially 
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occupied state, so that the number of unoccupied states is z^ — ^ — 1. We 
write a combined formula of the statistical weight for both the cases as 



Wi = ^77 (3) 
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where I = or 1 if rii/q is integer or noninteger, respectively; i — 1, 2, . . . , m. 

In the particular cases, q — 1 and q = rii, we have — = rii/q and I = 
so that Eq. (H) reduces to Fermi-Dirac and Bose-Einstein statistical weights, 
respectively, 

As one can see, the effective number of available single-particle states 
derived from Eq. 
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for fermions and bosons, respectively, is linear in rii. With the identification 
of the parameters, 1/q = A, and the redefinition, Z\ — > Zi — (1 — A), the 
statistical weight ([|) coincides with Haldane-Wu statistical weight ([[]), for 
the case of integer 7ii/q. Consequently, the obtained statistical weight (H) 
corresponds to a kind of fractional exclusion statistics. To verify whether 
(0) leads to Haldane-Wu distribution we obtain below the equation which 
governs statistical distribution. 



3 The distribution function 

Starting with Eq. (|3|), we follow usual technique of statistical mechanics to 
derive the associated most-probable distribution of n,. 

The thermodynamical probability is W = Y\Wi, and the entropy, S = 
kin W, can be calculated by using the approximation of big number of par- 
ticles, n! ~ n n e~ n for big n. Assuming conservation of the total number of 
particles, N = J2 n i and the total energy, E = J2 n i£i, variational study of S 
corresponding to an equilibrium state gives us 

5S = kJ2 
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5n t = 0, (6) 
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where a and /3 are Lagrange multipliers, and we have used 
I = for big rij. Using the notation k = 1/q and inserting a = —ii/kT 
and /3 = 1/kT (obtained via an identification of S, at q = 1, with the 
thermodynamical expression), we rewrite Eq. (^j) as 

(zi + (1 - - Kni) K £i - \i 11 

= exp - — , «; = 1, 7 
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To draw parallels with Haldane-Wu statistics below we make analytic 
continuation of the discrete parameter k assuming k G [0,1]. Under this con- 
dition, the derived distribution equation (0) does reproduce that of Haldane- 
Wu fractional exclusion statistics (Eq. (14) of ref. 0), with k = A. 

Below, we turn to consideration of properties and exact solutions of 
Eq. ©. 

In general, Eq. ([7|) can not be solved exactly with respect to rij. However, 
for k = 1 and k — > (m^ — > 1), it becomes linear in n« and gives Fermi and 
Bose distributions, respectively. Also, we note that for k = 1/2, 1/3, and 
1/4 the equation contains a polynomial of degree up to 4 so that it can be 
solved exactly for all these cases. 

A convenient expression for rij obeying Eq. ([?D is given by 

ni = , \- , (8) 

W(X) + K 

where we have redefined, rii/zi — > rii, x = exp[(£j — fj,)/kT], and the function 
w(x) satisfies 

(l + w) {1 - K) w K = x. (9) 

Remarkably, exclusons which are "close" to fermions can be described in 
terms of exclusons which are "close" to bosons. In fact, we note that Eq. ([?[) 
is invariant under a set of transformations, 

k — > 1 — /t, rii — > ~ n i, % — > —x, (10) 

for k 7^ 0,1. Therefore, if ) satisfies Eq. ((7|) then the function = 

— 7ij(— a;, 1 — «) satisfies the same equation. Thus, we obtain the following 
general relation 

Hi(— x, 1 — k) — — ni(x, k), k ^ 0,1. (11) 
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We see that the distribution rij of exclusons for, e.g., k = 1/200 ~ can be 
obtained from that of "dual" exclusons, with k = 1 - 1/200 = 199/200 ~ 1. 

The values k — 1 and k — > («n, — ► 1) are the only two points of 
degeneration of Eq. (|7|). Hence, any "deviation" from Fermi or Bose statistics 
is characterized by a sharp change of statistical properties, sending us to 
consideration of exclusons. Consequently, we can divide particles into three 
main types, genuine fermions, genuine bosons, and exclusons (k g]0, 1[), since 
their statistical distributions obey different non-degenerate equations. 

A fixed point of the map k ^ 1 — k is k = 1/2. Hence it represents a 
special case worth to be considered separately. In this case, Eq. (0) allows 
an exact solution and the result is (positive root) @ 

2 2 

= , = T7o- ( 12 ) 

(l + 4expfc^) 1/2 

This distribution represents statistics with up to two particles in state, q = 2 
(semions). 

We have obtained exact solutions (real roots) of Eq. (0) for k — 1/3 and 
2/3 which we write as 

* = /+A T i< (13) 

where 

V3 /3x\ 3 
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From Eqs. flED and (H) one can see how exclusons with k — 1/3 (upper 
sign) are related to exclusons with k = 2/3 (lower sign) that agrees with 
Eq. (p|) . Also, for k = 1/4 and 3/4 we have obtained the following exact 
solutions (positive real roots): 

4 . . 

m = = , (15) 

Pg-l/2 _ ~ + 3 ± gl/2 T 2 

where 
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z 2( z + 2 )] 1 /3 + [^ + 2 )2]V 3 )+l, z = a/3 — +1-1. (16) 
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Plots of nj(x) for various k are presented in Fig. 1, from which one can 
see that these exclusons behave similar to fermions. 
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Figure 1: Statistical distribution m as a function of x = exp[(&j — fi)/kT], 
for k — 1 (fermions), k = (bosons), k — 1/2 (semions, Eq. fll2|)), « = 1/3 
(Eq. (JTSj) , upper sign), k — 1/4 (Eq. (|15|), upper sign), k = 2/3 (Eq. 
lower sign), and k = 3/4 (Eq. (|T5|), lower sign). 



Distributions of exclusons can be obtained from a different approach, 
based on the canonical statistical sum which implies the mean number of 
particles, 

E Nx~ N 

n = ^l . (17) 

This formula gives (exact) Fermi and Bose distributions for q — 1 and g — > oo, 
respectively, while for arbitrary q > 1 the sum is 

x 1+q - (1 + q)x + q 
n = ~ (x^-l)(x-l) ' g=1/ " (18) 



In Fig. 2, we compare distributions ( |T8| ) with exact solutions shown in Fig. 1. 
One can see that deviations become considerable as k goes to smaller values. 
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However, we expect that near k = there should be a better correspondence 
since one approaches the other interpolation endpoint (bosons). We treat 



( p~8| ) as an approximate result which is useful since it gives a single simple 
distribution formula for all exclusons, k e [0, 1]. 

A connection between the two approaches requires a deeper study which 
can be made elsewhere. 
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Figure 2: Statistical distribution function of x = exp[(£j — /j)/ kT] , for 

various k; dashed lines represent the approximation to exact solutions 
(solid lines) shown in Fig. 1. 



4 Conclusions 

(i) The derived statistical weight and Haldane-Wu statistical weight 
([]]) lead to the same distribution equation (|7|); 

(ii) Haldane-Wu parameter A acquires a physical meaning of an inverse 
of the maximal allowed occupation number in state, X—l/q, similar to the 
inverse of the statistical factor as shown by Wu @; 
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(iii) Within fractional exclusion statistics, the generalized Pauli exclusion 
principle reads that a maximal allowed occupation number of identical par- 
ticles in state is an integer, q = 1,2,3,..., i.e. rii/zi < 1/A as formulated 
by Wu @. We stress that in our approach we use this principle as a basis 
to calculate statistical weight (§) rather than derive it a posteriori from the 
analysis of a statistical weight or distribution function; 

(iv) While Haldane-Wu parameter A is assumed to vary continuously, 
the statistical parameter k — 1/q runs over discrete set of values, k = 
1, 1/2, 1/3, . . . This may be an important difference since physically accept- 
able solutions of Eq. (^) may not exist for all values of k e]0, 1[, while 
k = 1,1/2,1/3,... guarantees a polynomial structure of Eq. (|7|), with physi- 
cally acceptable solutions; 

(v) The equation (0), which defines statistical distribution of exclusons, 
k g]0, 1[, has a remarkable symmetry ([II]) which allows to interconnect solu- 
tions Hi for k and 1 — k due to Eq. ([TT|). 

In summary, we have shown that Haldane-Wu fractional exclusion statis- 
tics finds a natural combinatorial and physical interpretation in accord to 
Eq. (D), and presented exact solutions of Eq. (0). 
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